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[i] We use measurements from the Lunar Exploration Neutron Detector (LEND) 
collimated sensors during more than one year of the mapping phase of NASA’s Lunar 
Reconnaissance Orbiter (LRO) mission to make estimates of the epithermal neutron flux 
within known large Permanently Shadowed Regions (PSRs). These are compared with the 
local neutron background measured outside PSRs in sunlit regions. Individual and 
collective analyses of PSR properties have been performed. Only three large PSRs, 

Shoemaker and Cabeus in the south and Rozhdestvensky U in the north, have been found 
to manifest significant neutron suppression. All other PSRs have much smaller 
suppression, only a few percent, if at all. Some even display an excess of neutron emission 
in comparison to the sunlit vicinity around them. Testing PSRs collectively, we have not 
found any average suppression for them. Only the group of 1 8 large PSRs, with area 
>200 km 2 , show a marginal effect of small average suppression, ~2%, with low statistical 
confidence. A ~2% suppression corresponds to ~125 ppm of hydrogen taking into 
account the global neutron suppression near the lunar poles and assuming a homogeneous 
H distribution in depth in the regolith. This means that all PSRs, except those in 
Shoemaker, Cabeus and Rozhdestvensky U craters, do not contain any significant amount 
of hydrogen in comparison with sunlit areas around them at the same latitude. 

Citation: Sanin, A. B„ et al. (2012), Testing lunar pennanently shadowed regions for water ice: LEND results from LRO, 
J. Geophvs. Res., 117, E00H26, doi:l 0.1 029/201 1JE003971. 


1. Introduction 

[ 2 ] More than 50 years ago, it was suggested that some 
areas near the lunar poles are sufficiently cold to trap and 
preserve for a very long time (^Gy) hydrogen bearing vola- 
tiles, either primordial or produced at the Moon via solar wind 
interactions or brought to the Moon as water ice by comets 
and meteoroids [ Watson et al., 1961; Arnold, 1979]. In the 
1990s, observations made by the bistatic radar instrument 
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onboard the Clementine spacecraft identified some anoma- 
lous regions on the Moon associated with permanently sha- 
dowed regions (PSRs) that were consistent with the presence 
of water ice [Nozette et al, 1996, 2001; Simpson and Tyler, 
1999]. These observations have been criticized as ambigu- 
ous [Fa et al., 2011]. Another attempt to test lunar PSRs for 
water was performed using neutron (LPNS) and gamma-ray 
(LPGRS) spectrometers onboard the Lunar Prospector mis- 
sion launched in 1998 [Feldman et al., 2004]. These instru- 
ments were able to create global elemental and neutron albedo 
maps and discover significant reductions of epithermal neu- 
tron fluxes in lunar polar regions poleward of 70° S and 70°N 
latitudes. This result was interpreted as an enhancement of 
hydrogen abundance in these areas [Feldman et al., 2001; 
Lawrence et al., 2006], 

[ 3 ] The advantage of nuclear methods over other remote 
sensing techniques is that they offer greater sensitivity to 
hydrogen (down to 100 ppm) plus the ability to probe to 
depths down to 1 m. Regional variations of neutron flux were 
observed by LPNS using its omni-directional sensors with 
poor spatial resolution relative to the smaller scales of most 
PSRs [Maurice et al., 2004]. 

[ 4 ] It has been shown that neutrons emitted from a planetary 
body with different energies are sensitive to soil composition. 
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The thermal neutron flux is dependent on the abundance of 
elements such as Fe, Ti, Gd, Sm. The fast neutron flux is 
proportional to the average atomic mass of the regolith and 
can serve as a marker that distinguishes between surface 
materials enriched with heavier elements (e.g., iron in the 
nearside of the mare basins) and lighter elements (aluminum 
in the far side highlands). These effects were clearly visible in 
global maps provided by LPNS [Feldman et al., 1998a, 
1998b; Elphic et al., 1998; Gasnault et al., 2001] and the 
recently processed LEND maps discussed by Lih’ak et al. 
[ 2012 ], 

[5] But the most intriguing result concerns deconvolution 
of regional variations of epithermal neutrons around lunar 
poles at a scale comparable with the size of large craters. It 
is thought that the distribution of epithermal neutrons in 
these regions is controlled by the hydrogen abundance. 
More than 10 years ago, LPNS found that the average 
neutron suppression in the epithermal energy range increases 
poleward from 70°, both in the South and the North, by 4-5% 
[Feldman et al., 1998b, 2001]. But the spatial resolution of 
these measurements was greater than most PSRs (FWHM = 
44 km [Maurice et al., 2004; Lawrence etal., 2006]), requiring 
model dependent data deconvolution of smaller surface fea- 
tures and the inference that hydrogen may concentrate at the 
bottom of permanently shadow craters [Eke et al., 2009, 
Elphic et al., 2007]. 

[6] From these results, the ensuing decade witnessed an 
ongoing series of sophisticated scientific investigations, dis- 
cussions, and speculations regarding the interpretation of 
these phenomena: was hydrogen localized in cold traps 
within PSRs and could hydrogen abundances be so high that 
it indicated the presence of water ice? Some publications 
indicated PSRs may contain up to 1.5% of water equivalent 
hydrogen presumably in the form of water ice from comets, 
while other research postulated that at least half of this con- 
tent could be hydrogen bearing volatiles created by solar 
wind interactions with the regolith, which subsequently 
migrated to the cold poles through micrometeoritic bom- 
bardment processes [Crider and Vondrak, 2000, 2003]. 

[ 7 ] Recent experimental observations from the Indian Space 
Research Organization (ISRO) Chandrayaan-1 mission, the 
NASA Lunar Reconnaissance Orbiter (LRO) and the Lunar 
Crater Observation and Sensing Satellite (LCROSS) missions 
add significantly to the discussion. The recent orbital radar 
observations of the Lunar North pole (MiniSAR onboard 
Chandrayaan-1) indicate differences in circular polarization 
ratios comparing the interior and exterior of some polar sha- 
dowed craters, which suggest the presence of water ice [Spudis 
et al, 2010]. Near- infrared spectral data have shown surficial 
hydrogen as either FLO or OFI is spread widely across the 
lunar poles, at depths up to ~ 1 mm (which is the approximate 
sampling depth of the near-IR instruments), and could have an 
abundance of ^100-1000 ppm FLO [Pieters et al, 2009; 
Sunshine et al., 2009; Clark, 2009], The LCROSS spacecraft 
observed the impact of its launch vehicle upper stage within 
the PSR of the south polar crater Cabeus and observed 5.6 ± 
2.9% of water by mass in the plume [Colaprete et al., 2010], 

[8] The most recent attempt using nuclear methods was 
made by the Lunar Exploration Neutron Detector (LEND) 
instrument onboard LRO [Mitrofanov et al., 2008, 2010a], 
LEND is multiband neutron spectrometer system similar to 
LPNS but equipped with a neutron collimator that significantly 


reduces incoming epithermal and thermal neutrons to the 
detector from non-nadir angles and yields an instrument foot- 
print ~10 km in diameter (Full Width at Half Maximum, 
FWHM, assuming the altitude of the spacecraft is about 
50 km). Furthermore, LEND’s field of view (FOV) is compa- 
rable with the size of many large PSRs of interest and provides 
additional capability to contrast neutron fluxes from the interior 
and exterior of shadowed regions. Preliminary analysis of 
LEND measurements was performed during the evaluation of 
potential impact targets for the LCROSS mission to quantify 
water content in the regolith. From these results it was shown 
that the physics of enhanced hydrogen regions was more 
complicated than was initially expected [Mitrofanov et al, 
2010b]. Some PSRs, like Cabeus and Shoemaker, are quite 
distinguishable with clear suppression of epithermal neutrons 
inside the crater compared with the fluxes in neighboring 
regions. These regional suppressions were interpreted as an 
indication of high hydrogen abundance and possibly water-ice 
deposits. But many other PSRs indicated no difference between 
the neutron fluxes in the PSR and neighboring sunlit areas. 

[ 9 ] In this paper, we present the first systematic analysis of 
neutron emission from PSRs. LEND data gathered during the 
primary part of the LRO mapping phase, which started on 
September 15, 2009, are used for this analysis and we have 
implemented several methods to compare fluxes in shadowed 
and sunlit regions at the same latitudes. A primary objective 
of this paper is to determine if all individual PSRs or a group 
of PSRs with small, medium or large areas contain a signifi- 
cant amount of water ice in the upper ~1 m of regolith in 
comparison with hydrogen content in the regolith of sunlit 
areas at the same latitude. The general smooth poleward 
increase of hydrogen content in the regolith observed by all 
neutron experiments at the Moon is not a subject of this paper. 

2. Instrumentation, Data Reduction and Methods 
2.1. Instrumentation 

[ 10 ] LEND is a neutron spectrometer with the ability to 
globally map neutron fluxes in different energy ranges 
including thermal, epithermal and fast neutrons. LEND also 
has the capability to measure epithermal neutrons with spatial 
resolution of about 10 km at the poles (see Mitrofanov et al. 
[2008, 2010a] for details). 

[ 11 ] The LEND mapping of regional variations of epi- 
thermal neutron flux across the lunar surface is supported 
using both an uncollimated detector Sensor for Epithermal 
Neutrons (SETN) with spatial resolution similar to LPNS 
(see global maps from this detector presented by Lih’ak et al. 
[2012] and the four Collimated Sensors for Epithermal 
Neutrons (CSETN 1—4). The CSETN 1-4 provide high spatial 
resolution observations, FWHM ~10 km [Mitrofanov et al, 
2008, 2010a]. The mapping of the SETN detector is similar 
to the measurement techniques applied to the LPNS instru- 
ment. The CSETN 1-4 3 He gas proportional counters are 
surrounded by a collimator that defines the aperture. The 
collimator sides and bottom are covered with cadmium and 
10 B [Mitrofanov et al, 2008, 2010a]. 

[ 12 ] Measurements made by the LEND collimated sensors 
are the primary object of this work, providing the unique 
possibility of measuring regions inside large craters with 
surface areas more than ^100 km 2 and to evaluate the local 
suppression or excess of epithermal neutron fluxes. 
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2.2. Data Reduction 

[ 13 ] In our analysis, we used LEND maps created from 
Planetary Data System (PDS) data including the period from 
September 15, 2009 up to December 15, 2010. A number of 
calibration processing procedures were applied to the raw 
LEND counting rates prior to mapping (see M. L. Litvak et al., 
LEND neutron data processing for the mapping of the Moon, 
submitted to Journal of Geophysical Research, 2012; and 
W. Boynton et al., High spatial resolution of epithermal neu- 
tron emission from the lunar poles: Constraints on hydrogen 
mobility, submitted to Journal of Geophysical Research, 
2012). This included: correction of counting rates for effi- 
ciency saturation effects, adjusting for long-term variation of 
Galactic Cosmic Rays (GCR), eliminating times of strong 
Solar Particle Events, separation of charged particle and neu- 
tron events, as well as temperature and altitude corrections. 

[ 14 ] Counting rates during orbital measurement in all LEND 
detectors are measured with time accumulation intervals of 1 s. 
For such short time intervals, the spacecraft flies over a surface 
distance of about ~ 1 .6 km. This is short enough to distinguish 
local landforms with sizes 10x10 km. For our analysis we 
have created polar map grids with a pixel size 1/2 degree 
longitude by 1/16 degree latitude. The latitudinal pixel size is 
comparable to the 1 .6 km distance. All counts measured by 
LEND collimated detectors were directly distributed into the 
corresponding pixels of this map grid with counting rate, 
exposure time and statistical error in each pixel. 

2.3. Testing Analysis of Permanently 
Shadowed Regions 

[ 15 ] Several analytical techniques are applied to the LEND 
high resolution neutron data in this study. The first is to 
create a detailed map of the distribution of epithermal neu- 
tron fluxes across polar regions and to search for local areas 
indicating significant Neutron Suppression Regions (NSRs) 
relative to background flux rates. (This analytical technique 
is described in Boynton et al. (submitted manuscript, 2012) 
and I. G. Mitrofanov et al. (Testing polar spots of water-rich 
permafrost on the Moon: LEND observations onboard LRO, 
submitted to Journal of Geophysical Research, 2012).) 
Localized NSRs can be compared with other available data 
describing properties of relief (such as elevation, roughness, 
and slopes directions), solar illumination and temperature 
distribution (day, night, average temperatures, model pre- 
dictions of possible ice depth). Another technique is to use 
known targets and landforms of the lunar surface and to 
correlate them with variations of epithermal neutron flux. 
For example, it is suggested that there may be a correlation 
between the orientation of crater slopes and LEND counting 
rates [see McClanahan et al, 2012], 

[ 16 ] In this paper we use the second technique and make 
estimates of the counting rates within known large PSRs 
(with area greater than or equal to the LEND footprint) and 
compare counting rates from the local neutron background 
measured outside PSR areas (sunlit regions). If our analysis 
shows that epithermal counting rates measured for PSRs are 
lower than local background, it would support the hypothe- 
sis that PSRs are major hydrogen rich areas on the lunar 
surface. Depending on the magnitude of this suppression, 
PSRs may even be considered as a primary target in the 
search for possible water ice reservoirs. 


[ 17 ] The LRO instrument payload includes the Lunar 
Orbiter Laser Altimeter (LOLA), which is specially designed 
for global mapping of lunar landscape and the creation of a 
new high resolution topographic map of the Moon [Smith 
et al., 2010]. After several months of operations the LOLA 
team was able to produce illumination maps of the lunar 
poles. The Moon’s spin axis is nearly perpendicular to the 
ecliptic plane, just ~1.5° away from the ecliptic pole. This 
creates natural conditions where sunlight never directly 
reaches the floor of some of the deepest craters in both North 
and South circumpolar regions. The new illumination maps 
provided by LOLA represent the percentage of time that the 
Sun was visible from a given point on the surface [Mazarico 
et al. , 20 1 1 ] . If the average solar illumination during the lunar 
precession period is equal to zero for some surface element, 
then it is defined as a part of a PSR. Clusters of contiguous 
pixels may be aggregated to define a single PSR. The full list 
of all PSRs both at the North and South poles consists of 
more than several thousand areas ranging from very small 
(<1 km 2 ) up to quite large regions of more than 1000 km 2 . 
For the South pole, the total number of PSRs is smaller than 
for the North, but on average the southern PSRs are larger 
with several tens of PSRs having areas more than 100 km 2 . 
The set selected for analysis in this study all have areas 
>100 km 2 . The list of tested PSRs is presented in Table 1. 

2.4. Local Background 

[is] To perform our analysis of PSRs we need to derive two 
key observational parameters: 1) the counting statistics across 
each PSR’s pixels and 2) the counting statistics acquired 
across neighboring sunlit areas, which we define as a local 
background. The latter value is based on some a priori infor- 
mation, and can be estimated by various methods. We assume 
the best option is to consider nearby sunlit areas having the 
same soil type and a large enough surface to get significantly 
large number of accumulated counts, much larger than ones 
derived from the PSR area. This selection criterion is neces- 
sary to ensure that possible variations between these two 
parameters are due to different H abundances and not due to 
regolith elemental variations. Our goal is to compare neutron 
flux inside PSRs with fluxes in surrounding sunlit areas to find 
local effect between shadowed and sunlit areas at the same 
latitude. The last condition is important because in earlier 
observations LPNS revealed the existence of extended neutron 
suppression which depends on the latitude. This extended 
neutron suppression shows a decrease in epithermal neutron 
flux as a function of increasing latitude poleward of 70° South 
and North. Therefore, the area selected for the estimation of 
local background should be located at the same latitude belt as 
the tested PSR. 

[ 19 ] Our first simple, conservative option for an estimation 
of local background around test PSRs is based on latitude 
analysis. We have averaged the LEND polar maps (see 
Figure 1) by longitude and created latitude band profiles of 
counting rates. All permanently shadow pixels have been 
excluded from averaging. The resulting profile shows a lati- 
tude dependence of counting rate for polar regions for sunlit 
areas only (see Figures 2 and 3). For each analyzed PSR 
we estimate the value of the local background using this 
profile, as the counting rate at the latitude belt containing 
the PSR itself. This option takes into account the latitude 
dependence of the neutron flux distribution across polar 


3 of 13 



Name of Crater-Contained PSR Area Exposure Min Lon Max Lon Min Lat Max Lat Local Local 

PSR or Nearest Crater (km 2 ) (sec) (deg) (deg) (deg) (deg) Suppression 1 Suppression 2 


E00H26 


SANIN ET AL.: TESTING LUNAR PSRs FOR WATER ICE BY LEND 


E00H26 


. O \© \0 \0 \© . © \© \0 \© \0 x© . © xO \© x® \© x© x© \0 \© x© x© \© x® \© \© x© \0 \© \© \© \© \0 \0 \© \© \© \0 \0 x© \0 \© \© \0 

Ax ©X ©N ©N ©N Ax ©x ©X o'- ©X ©X A\ ©x ©N ©X ©x ©X ©X ©X ©X ©X ©X ©X ©X ©x ©X ©X ©X ©X ©X ©x q\ ©x ©X ©X q\ ©x ©X AN qX ©X qX ©X ©X ©\ ©X 

^^c^^'Off 5 ^uS^t^vS^ , /Sfnvd , rj^r 4 vdvdTi-vdod^c^ONc 4 c^ON'Oo 6 ^c^ONvS , /Sr 4 ’-HOCNodfoovd<Nod 

_ xO xO xO x© _ xO X© X© X© X© _ X© xO xO X© X© X© X© X© xO xO X© "T 1 X© xO X© X© x© xO "n X© x© X© X© "Tl "Tl M ..© x© ~n “n li 'Ll 0 


x® ©N ©A ©A ©x n ® ©x ©X ©N ©X ©X 


xOxOxOxOxOxOxOxOx^Jn® 

©N ©N ©N ©N ©N ©N ©N ©N ©N ©N 


© ' O' O' O' O' X. O' O' O' O' O' O' O' O' O' O' O' X. O O' O' xO O' O' O' O' XO O' O' O' O' X O X © X O © ' O' X © X O \0 X O I 

^^^^qr^®^^oqqt v ;f0^tONa\^;'0(^o^^;o^'T;(NqTt;^i^GoiDfncNo x (^qo<^<A'o^o v ' 
'^fD^fDri‘fn(NAt^vo(N^ 6 nd'ri 9 ^dfn ir !diD-HDi pH . L'-H>DiD ,v ; (v ; ir !o 6 ro r i' : ^ a l , 0 | 

I Sill II 17 2 I v I 7171 - i<Noj- 222° 


X © -xO -x© x© X® x © X® X® X® X© x© x.© X® X© X® X® X© X© X® X® X© X® X® X© x® x® X® X® X® X® X® X® X® X® X® X® X® X® X® X® X® X® X® x® X® X® 

©X ©X ©X ©X ©X X ©X ©X ©X ©X ©X A\ ©X ©X ©X ©X ©X ©X ©X ©X ©X ©X ©X ©X o' O' o' o' o' O' o' ©X O' O' ©X O' o' O' O' o' o' O' o' o' o' o' 

N ^\q^vq^nt^ON^^^^D;Osoo'«q(>'tqoo(N , tfOD)-H^qh;(SONqqr^voD;ON"D'Doooocno^^tx;iD 

^^C'i^vd f ^uDvS^c^^ T juScosduS^cxisd^^'Ood^i>ONCxii^O'sv©od'^[^osuD<Ac4^©0'sodco©'Oc4od 

X® X® X® X® X® _ X® X® X® X® X® _ x® x® X® X® X® X® X® X® X® X® X® ^ X® X® ^ X® X® X® X® ^ X® X® X® x® ^ ^ ^ X® . ® -tl ^ . ® 

O'** ©^ ©^ ©^ O' ^x ® s ' ®^ ®'' ®^ ®^ Ax ® s ' ®x ©^ ©^ ©^ o' o' ©x ©x ©X ©X ® ©X ©X . © ©N ©x ©X ©X .© ©x ©X ©X ©X ® -v © \ © ©X o' \® x® \® \© ®x 

iri’TOcim^^iAh'too2'oo-ooRiAmmooocMo x ocoiA^ou'Oo'0'TRi'Ttf'o'o'-ntf'c'o'^M 


^ •t n a m (N , 


; ro <N SO CO SO ■ 

III I 


• O CO O CO . t 1 

2 V I I ’ 


<n ^ ^ ^ so 

— m (N so o — 1 


Trquooqqrn&oiririUF-^q^qwri^ooooriu^Ti^fNrHTt^ccqririfNjririqqiAriu^fN^'t 
hsdsduni/)^ 2 nri 2 Ttiricorj^sdoon^h^ddnci^iA(N^rirniA(Nsc)^(Arnri(Nri^--lnO'T 
00000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000 
I I I I I I I II I I I I I I I I I I I I I I I I I II I 


sq^sqsqoof^o — CNjt^^^^sqrqaN^t^OTt^oo^^o^t^sqossqcNjuDsqsqsqoqr^vT^-^sq-Trr-osaNr^ 
xoo'Aodrn^ri^A-^ri^Troo^rv^ATf-Tt-od^^d^Di-HiG^'so^fn'^risdiTiDjrnfNf^fnr-'drio^-T 
0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000^-00 
I I I I I I I II I I I I I I I I I I I I I I I I I II I 


0 "to\o^' 0 ^oinmooinMwvtiA\f'OfnOfninaiOooAm' 


mdSOOAASOOOOOO 


r)fS'tcn'noc-HL''-ioofnN'toom'Oso'si- , ooohfnt'h^h'DroiDmo\^owniDOO'H-H^r'M-H , to\fO 

so-'0\(Nasir)'ooMU'--i , tot^'^-'OM'^so r'lsoroco^-isoosinsotnsooo'sj- | m o ^ so a n vi 00 | 


OfT^OOf , )'TUM\t-in , t-'T(NOnOSCO(NinOmf f infnUA"tO\ONO'. SOahM'tCOOOONMOMsO 

r'r'i-^^oop'ri^^o 6 dfo^o 6 o^sdddfosdK^fnfndo 6 ^fi^ , TO\'t'>dt'rohr'fn i O'si ; o\ 0 \ 6 ^d 

n-hL''Ooo^iot--hoo w|i--orfiD^so'0 oinmn(NL' 0 ' 7 fso 7 isoL''ri n ri 0\ 1- so t-- i— 1 tj- 00 «/p | 


»hMoontvisoO'HmHMMh(Noooa\iAOmM^^'omwt , t'thoonvia\OrohOsfO 
oosqTtr^^ooor^rqos^oorqooooosqrqr^sqTt^oquDrqosr^'— ;oo'xt;osfNoooqos^^uD(Nrq 
i/irniorisddsdsdfnio^o 6 ddioro'oo 6 nririfno\dL'sdo 6 'troa''^d(N-- ' so 00 ro cd 00 00 r--^ 
sOOO-rOOOUUOO'T'AXTtTt^O'Tr-so-T'TOOOO'OO'-'^nCN^sO'tOh^sOOOMON'T 
ooomoo' 0 'ooso\'H'AooHOso\inoon^oh'HhfO'H'nnn'nfO'om-H>nnoooO'HpH(SMfri 


(N so so (N (N 

— Yi A 00 'Ll 


© 00 CO OO SO SO OS OS 

O N ' — 1 WH 


• in m cn ud m so m 


a-Hinmo»o\NOONNO\Mhoo' 0 (NOiaiso-HiD , TMO\ooiDsoa\omnmo(SOiD^o\o-Hn-HhooiD 

0 \Nr^iDc^' 0 (NrDqqqq^;'t-;fsq 7 t;(Noq o oiDiDqooqfDqfDqq'HrD(NrDq^ , tqq(Nqqoot^rD 

o\d | DdiDdDidiDdsdfn^'ci^iDd^f , idDiKt'o 6 Ri'-iodsd^ ; o\sdiDfnrndo 6 sd^ : fnri-Hodf<iridd 

h-HsoiDOO'oonM(NO»MT'Tnn-HO\o\couso'T^\fr f inciM(Nn( Nnnr - r -- H--rHooooo 

OOsouD^tf^mmromrori(N<NfN<NfN(N^^-— 1 ^— ir-H-— 1 — 


a> 

i ™ 8 s s 

1 1 2 1 % S 1 


Z w 

C- . 


“ w£ S 

■RSIS 


jj 0 A 
& o _ 

M CS u ” 


, -a to co <d tj 

! c c p 9 - o 




> -5 a £ 

030 ^>. 00 c 30 

iZwZU w _) ZZZ 


I 03 *2 O hT 

! o § a uh ffi ^ , 


& 

V) 

x - 

-m +3 ^ 00 1 

s S s a 

aJn 1 > U 
C3 O M 
— 3 a ) co co 

5 7 -2 £ , 

S > CU Clh , 


4 of 13 


‘Bold font indicates there are only three PSRs that have significant neutron suppression relative to local background. 
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Figure 1. LEND map of epithermal neutrons at the lunar poles. Raw counting rate maps are smoothed to 
decrease random fluctuations. Smoothing scale increases outward from the poles to produce similar count- 
ing statistics at all latitudes. Black spots and contours of the maps represent the boundaries of PSRs 
according to LOLA data [Mazarico et al., 2011], 


regions, significantly increases the counting statistics and 
decreases the statistical error in the estimation of local back- 
ground. However, it may not account for a possible longitude 
asymmetry of neutron flux suppression around the poles. 

[ 20 ] Another background estimation technique is more 
complex and uses a two-dimensional approach. In this case 
we again take the original polar maps of the LEND counting 
rate, but smooth with a Gaussian filter of pre-selected full 


width at half maximum (FWHM) much larger than size of 
PSRs, excluding in the smoothing process all known per- 
manently shadowed pixels. On this map all pixels within 
PSRs are filled by values predicted by the Gaussian 
smoothing process based on surrounding sunlit pixels. This 
approach provides an estimate of the local background 
counting rates from the sunlit vicinity around the selected 
PSR location. The disadvantage of this method is that we do 



Figure 2. Latitude profile of counting rate averaged over 1° latitudinal belts for North polar region. Sta- 
tistical uncertainties are shown as vertical error bars. 
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Figure 3. Latitude profile of counting rate averaged over 1° latitudinal belts for South polar region. Statis- 
tical uncertainties are shown as vertical error bars. 


not know by definition the best FWHM for the smoothing 
filter. Too small a filter size will provide an estimate of local 
background with significant uncertainty due to poor statis- 
tics. Too large a filter size will smooth out the effect of 
extended polar suppression of neutron flux. For testing large 
PSRs, we will use the smoothing filter with the distance 
scale of ~227 km = 7.5° of latitude (see below). The 
selected scale is the distance on the lunar surface from sub- 
LRO point to the horizon visible from the spacecraft at 
50 km altitude orbit (this distance is ^412 km) multiplied by 
a factor ^0.55. 

3. Testing of Individual PSRs 

[ 21 ] In this section we present results for testing epithermal 
neutron emission for large area PSRs with data from the 
LEND collimated detectors. We have taken time series data 
gathered from the four LEND collimated detectors for epi- 
thermal neutrons distributed across a map grid (as described 
in section 2) and estimated the average counting rate for 
each individual PSR having area larger than or comparable 
with the LEND spatial resolution. Taking into account that 
the LEND spatial resolution is a spot on the lunar surface 
with radius ^5 km, we selected only large PSRs with areas 
greater than 100 km 2 . There are 46 such PSRs, 29 in the 
South and 17 in the North polar regions. Their areas range 
from ~ 1 00 km 2 up to > 1 000 km 2 . Eighteen of them have area 
>200 km 2 . The two largest PSRs have areas ~1000 km 2 
(Shoemaker and Haworth). So, about half of the PSR set is 
significantly larger than the LEND spatial resolution. 

[ 22 ] All selected PSRs are presented in Table 1 together 
with their areas, coordinate box (minimum and maximum 
longitude/latitude limiting the PSR area) and the neutron flux 
suppression within PSR borders found using two methods of 
background estimation (see previous section). The values of 
suppression are presented as the difference, A, between the 
average counting rate estimated for the PSR and the counting 
rate measured for the surrounding sunlit areas where we 
assume ~100 ppm hydrogen in the regolith [Lawrence et al., 
2006]. In each case, we also estimated statistical errors to test 


statistical significance of measured differences. The local neu- 
tron suppression parameter is presented in percent according to 
the expression: 

P sup = ^xl00%, (1) 

where the value of 1 .7 is the averaged counting rate in the field 
of view of LEND collimated detectors outside of extended 
polar suppression region [see Lih’ak et al., 2012], 

[ 23 ] We consider the local suppression in a PSR to be 
significant if it differs from zero by more than three standard 
deviations or 3<r (this means that only 0.27% of these cases 
may be the result of a random fluctuation). For any PSR it is 
possible to find an upper limit (F) for its neutron suppression 
parameter, P sup . This limit is calculated to be the value of the 
suppression parameter required to make it significant at the 
3cr level. The upper limit for local neutron suppression is 
estimated from the expression: 

3 c 

r = ^ x 100%, (2) 


where 8 = y A p SR + A 2 hgd - is a total statistical error 

corresponding to the difference in counting rates between a 
PSR with uncertainty Aj>sr and local background with uncer- 
tainty A hgd ; the value of 1.7 is the same as in equation (1). 

[ 24 ] As indicated in Table 1 by bold font, there are only 
three PSRs that have significant neutron suppression relative 
to local background: 

[ 25 ] 1 . PSR in crater Shoemaker in South circumpolar area 
which has a suppression value in the range of —5.5% to 
—6.3% (depending on how neutron suppression is calculated, 
see above). These values correspond to 220-245 ppm of 
hydrogen, provided it is homogeneously distributed in the 
regolith (method of calculation is described in Mitrofanov et al. 
[2010b], 

[ 26 ] 2. PSR in crater Cabeus in South circumpolar area has 
suppression in the range— 14.4% to —14.9%. These values 
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Cross section along 50.0 longitude 



Figure 4. A cross-section in latitude at 50.0° longitude of Shoemaker crater for LEND (solid line), 
LOLA (dotted line) and Diviner (dashed line) data. The “LEND smoothing scale” corresponds to a 
FWHM of 10 km in the collimated detectors. Good correlation of minimum values in all three data sets 
is visible. Maximum hydrogen concentration in this crater is well correlated with crater typography and 
with minimum of observed average temperature of regolith. 


correspond to 460-470 ppm of uniformly distributed 
hydrogen. 

[ 27 ] 3. PSR in crater Rozhdestvensky U in North circum- 
polar area which has suppression in the range —10.6% to 
— 11.3%. These values correspond to 345-365 ppm of 
hydrogen concentration. This is the only PSR in the North 
circumpolar area with large suppression. Though the sig- 
nificance of suppression in this case is slightly less than 3cr 
we include it in this set because of its identification with a 
large crater in the North polar region. 

[ 28 ] None of the other tested PSRs have a statistically 
significant suppression of epithermal neutron flux relative to 
neutron flux from local background. From this result we 
may conclude that all polar PSRs do not appear to be pri- 
marily responsible for the extended suppression of epither- 
mal neutrons observed at the poles, and individual PSRs do 
not contain large amounts of water ice deposits to a depth of 
~1 m in the regolith. 

[ 29 ] The epithermal neutron flux suppression of the PSR in 
Cabeus crater is the largest statistically significant suppres- 
sion of neutron flux in the South circumpolar area and there- 
fore corresponds to the highest hydrogen abundance in the 
regolith. As discussed earlier, the PSR in crater Cabeus was 
the target for the LCROSS impact [see Colaprete et al. , 20 1 0] . 
From this experiment the estimate of water concentration in 
the soil at the impact site is larger (5.6 ± 2.9% by weight) than 
the value of hydrogen concentration found in this analysis 
(460^170 ppm). This difference may be explained, in part, by 
the fact that the LEND result represents an average over the 
total area of the PSR which is ~280 km 2 and much large than 
LCROSS impact site. Also, one may consider a double layer 
regolith model with a dry top layer a few tens of centimeters 
thick above an enhanced hydrogen lower layer. In this case, 


the observed suppression will correspond to increasing 
hydrogen concentration with increasing thickness of the 
overlying dry layer [Mitrofanov et al., 2010b]. 

[ 30 ] It is interesting to compare three large southern PSRs 
in craters Shoemaker (51 km diameter), Haworth (35 km 
diameter) and Faustini (39 km diameter). These craters have 
about the same floor depth compared to local terrain, their 
ages are about the same, and floor average temperatures are 
similar, less than 60 K for all three (data have been taken from 
the LOLA and Diviner PDS). But only PSR in Shoemaker 
manifests statistically significant neutron local suppression 
with possible concentration of 220-245 ppm of hydrogen in 
the regolith. To illustrate this result we present in Figure 4 a 
cross-section in the meridian direction of the crater relief, 
average surface temperature and measured neutron flux. Only 
in the case of Shoemaker crater does the maximum of 
hydrogen concentration correlate with both the crater floor 
and also with the minimum observed average temperature of 
the regolith. A slight shift of maximum hydrogen concen- 
tration from the crater’s center in the equatorial direction is in 
agreement with results obtained in McClanahan et al. [2012], 
Upper limits of neutron suppression, F, for the PSRs in 
Haworth and Faustini craters are equal to —4.3% and —7.9%, 
respectively. This means that in these two PSRs there is no 
hydrogen concentration larger than 195 ppm and 275 ppm, 
respectively. 

[ 31 ] In the case of the PSR in Cabeus crater, a good cor- 
relation also exists with the maximum hydrogen concentra- 
tion (according LEND measurements), minimum of average 
temperature (according Diviner data) and maximum floor 
depth (according LOLA data). This correlation is shown in 
Figure 5. It is a task for future work to understand why the 
PSRs in Shoemaker and Cabeus craters present strong 
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Cross section along -46.0 longitude 



Figure 5. A cross-section in latitude at —46.0° longitude of the Cabeus crater for the LEND (solid line), 
LOLA (dotted line) and Diviner (dashed line) data. “LEND smoothing scale” corresponds to a FWHM of 
10 km of the collimated detectors. The minimum of neutron flux (maximum of hydrogen) is well corre- 
lated with minimum of surface average temperature. Both minima are located at bottom part of the crater 
wall close to its floor. 


enhancement of hydrogen concentration, but PSRs in Haworth 
and Faustini craters do not. 

[ 32 ] Moreover, there are several PSRs with areas greater 
than 300 km 2 showing some epithermal neutron excess 
rather than suppression in comparison with local back- 
ground. There are two possible reasons for epithermal neu- 
tron enhancements: 

[ 33 ] 1. All found enhancements have statistical signifi- 
cance less than 3cr and, therefore, they may be a statistical 
fluctuations of measured signal; 

[ 34 ] 2. All tested PSR are located inside the general pole- 
ward neutron suppression with deepest suppression ~5%. 
This general suppression has been explained by general 
enhancement of H concentration in regolith with increasing 
of latitude. Therefore, some area above ±70° latitude may 
have higher epithermal neutron flux in comparison with 
surrounding regions at same latitude in case if this area is 
drier than its surrounding regions. In this case observed 
neutron enhancement must be less of equal to ^5% taking 
into account the statistical uncertainty of found value. 

[ 35 ] This leads us to the conclusion that large PSRs with 
sizes greater than the LEND FOV are not strong sources of 
neutron suppression of epithermal neutron flux as was pre- 
viously thought. Data show that large PSRs do not appear to 
be local spots of epithermal neutron flux suppression. In 
some particular cases they show significant but moderate 
depression of epithermal neutron flux in comparison with 
nearby values of neutron flux attributed to average extended 
polar suppression in sunlit areas. Assuming a homogeneous 
model for the regolith, we may estimate a higher hydrogen 
abundance of several hundred ppm, which does not require 


the presence of traces of dirty water ice. However, the depth 
distribution inside such PSRs is more complicated. 

4. Statistical Properties of Large PSRs 

[36] In this section we continue the discussion of possible 
hydrogen enhancement in PSRs, focusing not on the most 
significant cases discussed in the previous section but 
applying statistical methods to the whole ensemble of PSRs. 

[ 37 ] The total number of PSRs according to the latest 
LOLA illumination maps selected for this study is 344 for 
both lunar poles. There are 298 PSRs (summing north and 
south PSRs) having areas less than 100 km 2 , which is less 
than LEND spatial resolution (group #1), 28 PSRs with area 
between 100 and 200 km 2 , which are comparable with 
LEND spatial resolution (group #2), and 18 PSRs having an 
area larger than 200 km 2 , which are easily resolved by 
LEND (group #3). The minimum size of PSRs included in 
this data set is limited by the shortest available observation 
time made by LEND and corresponds to spots of a few 
kilometers in diameter. The minimum LEND data acquisi- 
tion time interval is 1 s. During this time LRO covers a 
distance of ~ 1.6 km on the Moon’s surface. 

[38] In this section we study and compare statistical 
properties of PSRs distributed into these three different 
groups paying special attention to groups 2 and 3 with large 
PSR surface areas. 

[ 39 ] In our first test we look at the distribution of the local 
suppression parameter (defined in expression 1) and its 
upper limit (defined in expression 2) as a function of expo- 
sure time (LEND observation time of a given PSR). For- 
mally speaking, the local suppression parameter (following 
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100 1000 10000 
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Figure 6. Local suppression/excess parameters (see expression (1)) measured for each large (area 
>100 km 2 ) southern PSR, presented as a function of exposure time accumulated during LEND observa- 
tions for each PSR. The dashed lines correspond to the 3a upper limits, T (see expression (2)), for excess 
(positive values) and suppression (negative values), respectively. 

from expression 1) may take both negative and positive level and higher, are located outside the dashed lines and 

values. In the first case, it indicates that the counting rate of values of lower significance are inside the lines. In these 

neutron flux measured inside the PSR is smaller than the figures one can see the largest PSRs (larger areas roughly 

average counting rate in the local sunlit vicinity. In the sec- correspond to longer exposure times) have high counting 

ond case, it shows that local vicinity emits fewer neutrons statistics and as a result small statistical errors. Only two 

than observed inside the PSR and local suppression pararn- cases of significant suppression are seen in the south 

eter in reality turns into local excess parameter. In Figures 6 (Shoemaker and Cabeus discussed in the previous section) 

and 7 we present distributions of the measured local sup- and one case approaching the 3<r significance level in the 

pression/excess parameters as a function of exposure time north (Rozhdestvensky U crater). For these PSRs the local 

for the selected south polar (Figure 6) and north polar suppression parameter ranges from —5 to —15%. Many 

(Figure 7) PSRs from groups 2 and 3 with areas >100 km 2 , other large PSRs, some with even less statistical error, have 

Local suppression parameter values, significant at the 3cr local suppression parameter within the same range, but well 



100 1000 
Exposure time, sec 


Figure 7. Local suppression/excess parameters (see expression (1)) measured for each large (area 
>100 km 2 ) northern PSR, presented as a function of exposure time accumulated during LEND observa- 
tions for each PSR. The dashed lines correspond to the 3a upper limits, T (see expression (2)), for excess 
(positive values) and suppression (negative values), respectively. 
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below the 3er significance level. Moreover, there are several 
large PSRs (with surface area 300-400 km 2 ) displaying 
excess neutron flux instead of suppression. This leads to the 
conclusion that ordinary randomly selected PSRs with area 
greater than 300-400 km 2 do not necessarily manifest local 
neutron suppression larger than a few percent in comparison 
with local surrounding sunlit areas. 

[ 40 ] In the previous section we have analyzed individual 
PSRs and were able to select only a few of them with sig- 
nificant local neutron suppression. On the other hand (see 
Table 1), there are many showing negative values of the 
local suppression parameter which could not be considered 
significant because their statistical uncertainty is too high. 
To decrease uncertainty one may combine PSRs in different 
groups and try to test the hypothesis that the distribution of 
counting rates within this given group is statistically differ- 
ent from the distribution of counting rates in the sunlit areas 
surrounding the PSRs included in the group. The counting 
rate in the sunlit areas is known very well (because of their 
large area and, as a result, enormous counting statistics) with 
the accuracy 10-100 times better than counting statistics for 
measurements inside the PSRs. The main point of the sta- 
tistical analysis suggested in this section is to reduce local 
background counting rate to the same uncertainty level as 
observed in the PSRs and then to compare it with distribu- 
tion of PSRs counting rates. The key value in such simula- 
tion is a PSR’s exposure time. It defines the statistical 
uncertainty showing how many counts have been collected 
by LEND detectors during observation of the given PSR. 
Using counting rate (counts per seconds) in the sunlit area 
and PSR exposure time T, (seconds) one may simulate the 
counting rate S, (counts per seconds) inside a PSR using 
Poisson statistic by the following formula: 

S, = Cf m /T u (3) 

Where Cf m is the randomly selected number of counts from 
the Poisson distribution with mean value B t x T r 

[ 41 ] From expression (1) we may estimate a simulated 
local neutron suppression parameter as 

Pf m = x 100%. (4) 

Following this method it is possible to estimate the local 
neutron suppression parameter for each PSR from the local 
background counting rate which is weighted with the expo- 
sure time of the PSR (thereby weighted with the real statistical 
uncertainty attributed to the PSR) and randomized using 
Poisson statistics. 

[ 42 ] The next step in our approach is to ensemble PSR 
groups for such analysis. We limited our selection to groups 
1-3 already introduced in the beginning of the section. This 
selection criterion combines PSRs according to two impor- 
tant properties: 1) corresponding PSR surface area and 
LEND field of view, 2) combining PSRs with similar sta- 
tistical uncertainties. 

[ 43 ] For each group of PSRs we have created two distribu- 
tions. The first represents the measured local neutron suppres- 
sion parameter (presented in Table 1). The second one is a 
Monte Carlo simulation of the local neutron suppression 
parameter using equations (3) and (4). The comparison between 


these two distributions may indicate if we are able to explain 
the observed diversity of PSR suppressions (within given 
groups 1-3) just by statistical variations of counting rates 
measured at the nearby sunlit areas. To make this comparison 
robust we have used two well-known statistical criteria. First 
is the Student’s t-test [e.g., Press et al., 2007]. It tests the null 
hypothesis that the means of two normally distributed 
populations are equal. Here, we suppose that our two tested 
samples of real and simulated PSR suppressions belong to a 
normal distribution, but may have different mean values. 
Using sample parameters such as mean value and variance 
we may create a t-statistic parameter and check the proba- 
bility that it follows a Student’s t distribution. Another cri- 
terion is based on the Kolmogorov-Smimov test or K-S test 
[e.g., see Press etal, 2007]. It is a nonparametric test that can 
be used to compare two samples by quantifying the distance 
between the empirical distribution functions of two samples. 
It tests the null hypothesis (calculates the probability) that the 
samples are drawn from the same distribution. 

[ 44 ] The real and simulated distributions of local neutron 
suppression parameter for the group #1 with 298 smallest 
PSRs are presented in Figure 8. Even from a visual com- 
parison, it is clear that two distributions are very close to 
each other. This is confirmed by both test criteria. The K-S 
test shows that measured and simulated distributions for 
group #1 display the same statistics with the probability of 
97%. Verification based on the Student’s t-test shows that 
these distributions have the same mean value with a proba- 
bility greater than 60%. So, we may conclude that statistical 
criteria cannot distinguish between neutron emissions inside 
PSRs with area <100 km 2 and in the sunlit areas around 
them. 

[ 45 ] Group #2 includes 28 intermediate PSRs with surface 
areas between 100 km 2 and 200 km 2 with the average area of 
~132 km 2 . These PSRs are comparable with the LEND 
footprint and could be resolved by LEND provided they 
differ from the local vicinity and have enough counting 
statistics to measure this difference. The comparison 
between samples of measured and simulated counting rates 
for this group (presented in Figure 9) shows that the proba- 
bility of coincidence is lower than was found for group #1, 
but is still too large to conclude that these distributions are 
statistically different. The value of the t-statistic parameter 
estimated for these samples may be randomly selected from 
the Student’s t-distribution with a probability of about 30%, 
which is enough to make the statement that these distribu- 
tions could have the same mean values of local suppression. 
The K-S test shows that shapes of the two distributions can 
be fitted to each other with a probability of 20%, which is 
also too high to distinguish them with high significance. We 
may only summarize that medium-sized PSRs from group 
#2 tend to emit fewer neutrons then surrounding sunlit areas, 
but this difference is small in comparison with statistical 
uncertainty and may be explained by random process of 
counting statistic, but not a physical difference in the 
mechanism of neutron emission. 

[46] Finally, group #3 includes the 18 largest PSRs with 
individual surface area greater than 200 km 2 and with 
average area ~407 km 2 , which is significantly larger than 
the area of LEND footprint. The comparison between 
observed and simulated distributions of counting rates for 
group #3 is presented in Figure 10. It is seen that their shapes 
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Figure 8. Distribution of counting rates measured inside small PSRs (surface area <100 km 1 ) shown in 
solid black. The distribution of simulated counting rates for these PSRs is shown by the dashed line. 


and mean values are quite different. If one uses the K-S test 
to check if these samples are taken from the same distribu- 
tion the most probable answer will be “no.” The probability 
to coincide is only about 1% and that is very close to the 30- 
level (probability = 0.27%) commonly used to define an 
unlikely event. The Student’s t-test also shows that the 
probability that these samples belong to a normal distribu- 
tion with same mean value is fairly low, ~0.7%. Analyzing 


these results we may state that the neutron flux measured in 
large area PSRs is probably less than the neutron flux in 
nearby sunlit areas. In our statistical analysis we also esti- 
mate the average local suppression, which would be neces- 
sary to achieve consistency between measured and simulated 
distributions. We shift the simulated distribution of the local 
neutron suppression parameter for the large PSRs until both 
statistical criteria show maximal probability of coincidence 



Figure 9. Distribution of counting rates measured inside medium-sized PSRs (with surface areas 100 — 
200 km 2 ) shown by solid black line. The distribution of simulated counting rates for these PSRs based 
on hypothesis about zero local suppression is shown by the dashed line. 


11 of 13 




E00H26 


SANIN ET AL.: TESTING LUNAR PSRs FOR WATER ICE BY LEND 


E00H26 


0.14 

0.12 

0.10 

2 - 

CO 

c 

CD 

0.08 

n 

co 

o 0.06 
Q. 


0.04 


0.02 


0.00 


Figure 10. Distribution of counting rates measured inside large PSRs (surface area >200 km 2 ) shown by 
solid black line. The distribution of simulated counting rates for these PSRs based on assumption of zero 
local suppression is shown by the dashed line. The distribution of simulated counting rates for these PSRs 
based on assumption of best fit local suppression is shown by the dotted line. 
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between the shifted distribution and the distribution of real 
measured suppression. The best shift value corresponds to a 
local suppression of —2.2 ± 1.5%. The result of this analysis 
is shown in Figure 10. Summarizing, one may say that large 
area PSRs belonging to group #3 display local suppression 
and may be treated as spots with enhanced hydrogen con- 
tent, but the value of the average local suppression is too 
small to make definite statement that large PSRs are very 
well localized and isolated spots with very high hydrogen 
content which may interpreted as water ice. 

5. Conclusions 

[ 47 ] Both analyses of individual PSRs and studies of groups 
of PSRs have shown that these spots of extreme cold at the 
lunar poles are not associated with a strong effect of suppres- 
sion of epithermal neutron emission. To estimate the local 
effect of suppression in PSRs, we have compared the count 
rates to the sunlit vicinity around them. We found only three 
large PSRs, Shoemaker and Cabeus in the south and 
Rozhdestvensky U in the north, which manifest significant 
neutron suppression, about —5.5% to —14.9%. All other 
PSRs have much smaller suppression, no more than few per- 
cent, if at all. Some PSRs even display a positive deviation, 
implying an excess of neutron emission in respect to the sunlit 
vicinity around them. 

[48] Testing PSRs collectively, we have not found any 
statistically significant average suppression. Only the group 
of 18 large PSRs, with area >200 km 2 , shows a marginal 
effect of small average suppression of about 2%, with small 
statistical confidence. LEND will collect more data and 
provide better counting statistics for PSRs, but even now the 
data is enough for definite conclusion that PSRs at both 
poles are not reservoirs of large deposits of water ice, which 


would be seen as depressions in the emission map of lunar 
neutrons. The major effect of extended suppression of neu- 
trons at the lunar poles is associated with sunlit surfaces, and 
one has to test the local neutron suppression regions (NSRs) 
outside the permanent darkness of PSRs. 

[ 49 ] On the other hand, the exceptional cases of Shoemaker, 
Cabeus and Rozhdestvensky U show that large enhance- 
ments of hydrogen exist in some particular PSRs. Compari- 
son of data for neutrons, surface altimetry and surface 
temperature shows very good agreement between the shadow 
in Shoemaker and emission of epithermal neutrons (see 
Figure 4). The presence of water ice has also been experi- 
mentally confirmed for regolith in Cabeus (5.6 ± 2.9%) by 
direct measurements of the plume from the LCROSS impact. 

[ 50 ] The physical processes for the trapping of volatiles is 
not as straightforward as previously thought, with the trap- 
ping of water molecules from the exosphere in cold spots of 
permanent shadow. One must study more complex physical 
models, where all PSRs are not the principal depository of 
lunar water and where there should be some physical dif- 
ference between PSRs, which make some better candidates 
for higher water content compared with the regolith in 
nearby sunlit regions. 

[51] Acknowledgments. We wish to thank the International Space 
Science Institute (ISSI, Bern, Switzerland) for the support of research 
(included in the framework of international team “Nuclear Planetology” in 
2007-2010) presented in this paper. 
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